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SYSTEM AND METHOD FOR A LITHOTRIPTER 

CROSS-REFERENCE TO RELATED APPLICATION 

This application claims priority to co-pending German Patent Application No. 102 34 
5 144.3, which was filed on July 26, 2002 and is incorporated herein by reference. 

FIELD OF THE INVENTION 

The invention generally relates to the application of ultrasonic waves to a target object in 
a body, and more particularly, to a system and method for a lithotripter that can effectively apply 
1 0 such ultrasonic waves. 

BACKGROUND OF THE INVENTION 

Lithotripters are nowadays widespread as medical devices for destroying concrements, for 
example kidney stones, in the body of a patient with the help of focused Shockwaves. Such 

15 lithotripters are commercially offered by the Applicant, for example under the designation 
"Domier Lithotripter S" or "Domier Compact Delta In all of these devices, the concrement 
must be located prior to the beginning of the treatment so that the patient can be positioned with 
the help of a movable stretcher such that, for example, his kidney stone is within the focus of the 
Shockwaves that are generated with the help of the Shockwave generator of the lithotripter. This 

20 initial "adjustment," i.e. positioning of the patient, is normally carried out with an imaging 
locating device, for example an imaging ultrasonic scanner or an X-ray locating device. In 
addition to the initial locating of the concrement for the first time before the beginning of the 
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ESWL (extracorporeal shockwave lithotripsy) treatment or ESWT (extracorporeal shockwave 
therapy) treatment, this device also serves to continuously monitor the position of the concrement 
during the treatment in order to make sure that it has not slipped within the patient's body or 
traveled to another place or that the patient has not moved on his stretcher, so that the 
5 concrement might no longer be in the focus of the Shockwaves. For a comprehensive survey of 
technical and medical aspects of ESWT and of the devices used in lithotripters, reference can be 
made to the book "ESWT and Ultrasound Imaging of the Musculoskeletal System," Steinkopff- 
Verlag, Darmstadt, 2001, ISBN 3-7985-1252-3. 

As a general rule, the imaging locating device is also used to monitor the success of 

10 lithotripsy in the course of the treatment, i.e. fragmentation of the target object. Since a treatment 
typically takes about 30 minutes, it is usually not feasible to complete the treatment in a 
continuous manner, because of the excessive radiation dosage of X-rays. Therefore, a treatment 
is typically conducted during intervals of 3 to 5 minutes. If the target object shifts during the 
treatment intervals, for example due to a movement of the patient, the patient's body will be 

15 loaded by Shockwaves, until the next control inspection, without the target object being further 
fragmented, since it is no longer within the shockwave focus. 

In lithotripters in which an ultrasonic scanner is used as the imaging locating device, the 
scanner may be used continuously for visualizing the target object, but the locating operation is 
often much more difficult than in radiograms, among other things, because of the image quality, 

20 so that even experienced medical personnel will often have difficulties in recognizing the target 
object or even in evaluating the degree of fragmentation of the object. Moreover, the use of such 
imaging locating devices not only for the "initial adjustment" of the patient, but also for "hit 
control" in the course of the treatment has a fundamental drawback. That is, although typically 
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the position of the target object relative to the Shockwave focus can thereby be controlled, for 
example by a purely geometrical quantity, the effect of the Shockwaves on the target object itself 
cannot be detected. Fundamental problems such as an inadequate coupling of the shockwave 
device to the patient's body, a vignetting of Shockwaves (e.g., by ribs), etc., might therefore be 
5 detected very late, such as when no effects become visible on the target object in the course of 
the treatment. 

To address the foregoing problems, special ultrasonic methods, particularly ultrasonic 
Doppler methods, for continuous hit control have been suggested in existing methods, based on 
the assumption that, when hit, the concrement in the human body will perform a macroscopic 

10 movement due to pulse transmission from the shockwave. When the target object is exposed to 
ultrasonic waves and when the ultrasonic waves reflected on the object are measured, this; 
macroscopic movement will be expressed in a Doppler shift of the frequency of the reflected; 
waves. Lithotripters equipped in this way, where the evaluation of the received ultrasonic waves 
includes a Doppler analysis, are known, for example, from EP 0 367 1 16 Bl, EP 0 548 048 Bl, 

15 and DE 44 46 192 Al. These devices have in common that the ultrasonic transducer emits 
ultrasonic waves into the body and that the ultrasonic waves reflected by the body back to the 
ultrasonic transducer are sensed, with a Doppler signal unit, which consists of the emitted and 
received ultrasonic waves, generating and evaluating a Doppler signal, wherein essentially the 
amount of a frequency shift of the reflected signal relative to the emitted waves is calculated and 

20 the hit accuracy is inferred therefrom. 

However, the foregoing approach has various drawbacks. For example, it is known that 
the Doppler signal includes artifacts, especially near its zero point in time (i.e., directly after 
emission of the shockwave). To prevent the measurement of such artifacts, the lithotripters 
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according to EP 0 367 1 16 Bl and EP 0 548 048 Bl are, for example, equipped with means for 
the time synchronization between the Shockwave generator and the Doppler signal unit, which 
makes these devices complicated and expensive. Moreover, in the vicinity of strong scatterers, 
e.g. a concrement, there may be the so-called mirror artifact in the Doppler spectrum. Due to 
5 double reflection, the movement of a scatterer is additionally recorded in the opposite direction. 
This leads to an additional amount, which corresponds to the useful signal reflected at the zero 
line. Due to the briefness of the process and the strong artifacts, velocity-time curves, and thus 
hit information, cannot easily be assigned to the spectra. 

Therefore, there is a need in the art to provide a system and method for a lithotripter that 
10 determines and displays information about hit control and disintegration control on the basis of 
received ultrasonic waves, without the need for performing a Doppler analysis, such as that 
described above, which requires complicated apparatus and is prone to errors. 

SUMMARY OF THE INVENTION 

15 The present invention is generally directed to a system and method for a lithotripter. In 

one aspect, the invention can include a system for a lithotripter, which can include a Shockwave 
generator and an ultrasonic transmitting/receiving unit that communicate with an evaluating unit. 
The ultrasonic transmitting/receiving unit can include an ultrasonic transducer that can emit 
ultrasonic pulses and receive ultrasonic waves reflected from a target area of the Shockwave 

20 generator. Furthermore, the evaluating unit can determine a correlation coefficient of the time 
correlation between the ultrasonic waves and the ultrasonic pulses and also provide a signal 
related to the correlation coefficient. 
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These and other aspects of the invention will be described further in the detailed 
description below in connection with the drawings and the claims. 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 The accompanying drawings are incorporated into and form a part of the specification for 

the purpose of explaining the principles of the invention. The drawings are not to be construed 
as limiting the invention to only the illustrated and described examples of how the invention can 
be made and used. Further features and advantages will become apparent from the following, 
and more particular description of the invention as illustrated in the accompanying drawings, 
10 wherein: 

Fig. 1 is a simplified illustration of major components of a lithotripter according to. 
exemplary embodiments of the invention. 

Fig. 2 is a schematic view of the evaluating unit and its connections in the lithotripter 
according to exemplary embodiments of the invention. 
15 Fig. 3 shows an example of two successively received ultrasonic wave echoes that were 

reflected on a static system following the miss of a Shockwave. 

Fig. 4 shows an example of two successively received ultrasonic Shockwave echoes that 
were reflected on a system in motion following a hit of the Shockwave. 

Fig. 5 shows a schematic example of the variation of the correlation coefficient with time 
20 upon a hit of a Shockwave. 

Fig. 6 shows a schematic flow diagram for illustrating the major sections of the method 
according to exemplary embodiments of the invention. 
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Fig. 7 shows a schematic flow diagram for illustrating method steps of the first main 
section illustrated in Fig. 6. 

Fig. 8 shows a schematic flow diagram for illustrating method steps of the second main 
section according to Fig. 6. 
5 Fig. 9 shows a schematic flow diagram for illustrating method steps of the third main 

section according to Fig. 6. 

Fig. 10 shows a schematic example of a relaxation-time variation control during a 
shockwave therapy. 

1 0 DETAILED DESCRIPTION OF THE INVENTION 

The illustrative embodiments of the present invention will be described hereinafter with 
reference to the drawings, wherein like elements and structures are indicated by like reference 
numbers. According to exemplary embodiments of the invention, an ultrasonic transducer of a 
generic lithotripter is designed to emit pulsed ultrasonic waves, and an evaluating unit is 

15 designed to determine a coefficient of the time correlation between reflected ultrasonic waves, 
which are assigned to successively emitted ultrasonic pulses, and to provide an associated 
correlation coefficient signal. This inventive design of the lithotripter exploits the fact that the 
more static the reflected system is, for example the body portion in which the target area of the 
shockwave generator is positioned and in which the ultrasonic waves that are emitted in pulses 

20 are reflected, the greater the correlation between successively received "echoes" (i.e., reflected 
ultrasonic waves). 

The following observation can be illustrative in this regard. When the reflecting system, 
for instance a kidney stone in a kidney of a patient to be treated, has not changed from one 
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emitted ultrasonic pulse to the next one, for example, for the reason that no Shockwave has been 
emitted by the Shockwave generator or for the reason that although such a Shockwave has been 
emitted, it has resulted in a "miss" (i.e. the kidney stone has been missed), the first echo (i.e. the 
reflection of the first emitted ultrasonic pulse) and the second echo (i.e. the reflection of the 

th 

5 second emitted ultrasonic pulse) are largely identical. When the time curve of the i received 
ultrasonic echo is generally designated as e\(t) and thus the time curve of a subsequently received 
ultrasonic echo, which represents the reflection of the (i+k) th emitted pulse, as ei+k(t), the time 
correlation between the two reflected ultrasonic waves can be defined as: 
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K i>k : 


T 

fei(t)e i+k (t)dt 
t, 


1/2 


1/2 


Je?(t)dt 


Je? +k (t)dt 


Eq. 1 


In this regard, the lower integration limit T\ and the upper integration limit T2 define the time 
window in which the received echo signals are evaluated. Ti and T2 are here chosen such that 
the echo originates from the target area of the lithotripter. The mean value of the time window, 
15 (Ti + T 2 )/2, corresponds to the distance between the ultrasonic transducer and the target area of 
the lithotripter. On the other hand, (T2 - T\)/2 defines the size of the volume used for the 
evaluation. 

In the above formula of Eq. 1, the denominator represents a standardization of the 
correlation coefficient Kj k . In the first considered case of a miss, where the reflecting body area 
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of the patient remains largely static, the ultrasonic echoes ej(t) and ej+k(t) are identical. Hence, 
their correlation Kj 5 k is substantially 1 . However, due to noise that usually can not be excluded 
entirely in practice, ej(t) and ei + k(t) may not be exactly identical, so that K ijk < 1 is typically 
applicable in practice. 

5 This situation will change entirely if instead of a miss a hit is achieved with a focused 

Shockwave, for this leads to different effects, each acting on the time curve of the ultrasonic 
waves received. An interaction mechanism between the Shockwave and the concrement leads to 
a macroscopic movement of the concrement, especially when it has already been fragmented in 
part or even for the most part. Especially at the beginning of the treatment, when the concrement 

10 is still unfragmented to a substantial extent, other effects are prevailing, namely above all the 
shooting of fragments out of the concrement and the occurrence of cavitation bubbles around the. 
concrement in the case of a successful hit. All of these mechanisms lead to a dynamic movement 
of areas in the patient's body on which the emitted ultrasonic pulses are reflected and have the 
overall effect that the stronger the interaction between the Shockwave and the concrement and the 

15 more the concrement has already been fragmented, the smaller the correlation becomes between 
ultrasonic echoes corresponding to the reflection of successively emitted ultrasonic pulses. 
Hence, the correlation coefficient signal provided by the evaluating unit includes the desired 
information on the strength of the interaction between the emitted shockwave and the target 
object. 

20 As can be gathered from the terminology used above and regarding the two ultrasonic 

echoes ej(t) and ej+k(t), two received ultrasonic waves can be used for determining the time 
correlation coefficients, which belong to emitted ultrasonic pulses that were not emitted directly 
one after the other, but are separated by further emitted ultrasonic pulses. However, it is 

8 


• ' Atty. Docket Na 08340.Lithotripter 

(Client Docket No. PA34859US-03476/wl) 

preferable that the evaluating unit be designed to determine the correlation coefficient with the 
help of reflected ultrasonic waves that are directly assigned to successive emitted ultrasonic 
pulses. Hence, in the above terminology, this corresponds to the case k - 1 and permits the 
maximal exploitation of the information received by the evaluating unit. In this , instance, a look, 
5 for example,, at the first three ultrasonic pulses that have been emitted after the generation of a 
Shockwave will reveal that the three ultrasonic pulses are reflected in the target area of the 
Shockwave generator, received as echoes 1 to 3 by the ultrasonic transducer, and evaluated in the 
evaluation unit, which determines the correlation between the first and the second echo, then 
between the second and the third echo, etc. Following each of the correlation coefficient 

10 calculations, an associated correlation coefficient signal is output by the evaluating unit. 

As for the further processing of the hit information contained in the correlation coefficient - 
signal, numerous procedures are possible. In one design of the lithotripter according to 
exemplary embodiments of the invention, the lithotripter further comprises a display device 
connected to the evaluating unit for displaying the time curve of the correlation coefficient. The 

15 correlation coefficients are represented in a coordinate system, where the abscissa is the time 
since the emission of the last Shockwave and the ordinate is the determined correlation 
coefficient. In the case of a miss, the time curve of the correlation coefficient has a constant 
value of about 1 over almost the entire period of time between two successive Shockwaves. By 
contrast, in the case of a hit, a "sag" in the time curve of the correlation coefficient is observed, 

20 due to the described effects in the target object that influence the reflected ultrasonic echoes. 
According to general experience, after about 50 to 100 ms, the fragments and the concrement 
itself will settle down again, and cavitation bubbles will have disappeared again after that period, 
so that successive ultrasonic echoes will again resemble one another more and more. The 


Atty. Docket No. 08340.Lithotripter 
(Client Docket No. PA34859US-03476/wl) 

correlation coefficient will then relax towards the value 1 . Thus, the medical personnel operating 
the lithotripter can assess whether there is a miss or a hit based on the display, according to 
exemplary embodiments of the invention. 

In addition, or as an alternative, the evaluating unit may be designed to output an error 
5 signal if, after the emission of a Shockwave, the minimum value of the correlation coefficient 
does not fall below a predetermined first threshold value. It follows from the above explanation 
of the impacts of a hit on the correlation coefficient that the time curve thereof shows a sag that 
becomes more pronounced as the hit becomes better (i.e., the stronger the interaction was 
between the Shockwave and the target object). Accordingly, the sag in the correlation coefficient 

10 in the case of a precise hit should be more pronounced. In other words, a miss or a "partial hit" 
can be recognized due to the fact that in the time-resolved display of the correlation coefficient, 
which has been described above, the minimum value of the coefficient does not fall below the 
first threshold value. The output of an error signal in this instance makes it possible to further 
process the miss information automatically, without an intervention of the medical personnel 

1 5 being needed. 

As an another addition or alternative, the evaluating unit may be designed to output an 
error signal if, after emission of a Shockwave, a relaxation time of the correlation coefficient falls 
short of a predetermined second threshold value, since a hit is also expressed by the feature that 
the system, consisting of concrement, fragments, and cavitation bubbles, needs a specific 
20 minimum time to settle down again according to a minimum relaxation time of the correlation 
coefficient (i.e., the width of the sag when the correlation coefficient is plotted in time). More 
specifically, in the course of the whole treatment, the relaxation time can be represented in the 
form of a continuous display to trace the development of the relaxation time. However, with a 
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successful continuing destruction of the target object, the relaxation time of the correlation 
coefficient should increase continuously. By contrast, if the relaxation time decreases, this is a 
strong sign that the concrement has moved out of the shockwave focus, either due to a movement 
of the concrement (e.g. the movement of a kidney stone within the kidney of the patient) or due 
5 to a displacement of the patient himself. If desired, the screen may also be used as a display 
device on which the images supplied by the imaging locating device are shown. 

The relaxation time can be determined in different ways. For instance, the evaluation unit 
can be designed to determine the relaxation time by adapting a Gauss curve to the time curve of 
the correlation coefficient. Depending on the curve shape occurring in practice, it may also be 

(— t/T 

10 expedient to adapt a function of the form 1-A V r; to the time curve of the correlation 
coefficient, wherein Tr is the relaxation time to be determined by adaptation. Likewise, it is 
possible in some embodiments of the invention to define the relaxation time T R via threshold 
values, i.e. as the difference of those points of time at which the time curve of the correlation 
coefficient falls below a specific threshold value for the first time or, in the relaxation phase, will 

15 then rise again above the threshold value. Such curve adaptations ("fitting") and threshold value 
analyses can nowadays be carried out by numerous data processing systems rapidly and reliably, 
so that a real time control of the shockwave therapy is possible with the help of the evaluating 
unit. 

In an advantageous development of some exemplary embodiments of the invention, a 
20 lithotripter may comprise an alarm device connected to the evaluating unit, which is fed with the 
error signal, and in this case, the alarm device can be expediently designed to output an optical 
and/or acoustic alarm. Thus, the alarm device can, for example, emit a warning sound or a 
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warning light as soon as it follows from the supplied error signal that the minimum value of the 
correlation coefficient does not fall below the first threshold value or that the relaxation time 
thereof falls below the second threshold value, which indicates a miss in each case. In addition, 
or as an alternative, the alarm may be triggered when the relaxation times measured show a 
5 negative trend, since in the case of a disintegration, a rising relaxation time must be assumed. 

Alternatively, or in addition, the Shockwave generator may be connected to the evaluating 
unit and designed to stop or continue the generation of Shockwaves as a function of the error 
signal. An automatic switching off of the Shockwave generator is thereby accomplished, i.e. 
when the target object is apparently no longer within the Shockwave focus. Hence, 

10 independently of the attention or of the reaction time of the medical personnel operating the 
lithotripter, an unnecessary loading of the patient's body by Shockwave misses can be avoided. 

The above-mentioned threshold values, which, when exceeded or fallen below, have the 
effect that the evaluating unit outputs error signals, can be preset in the evaluating unit as a rule 
and both the technical data of the Shockwave generator (for example, its performance) and typical 

15 data of the patient can be taken into account in this process. Expediently, however, the 
evaluating unit can comprise adjusting means for adjusting the first and/or second threshold 
value. In this instance, the threshold values can also be set individually at the beginning of the 
treatment, for example, in consideration of the arising signal noise. The threshold values can 
also be determined automatically on the basis of the signal itself, since among the signals 

20 acquired prior to the beginning of the shockwave therapy, it is possible to determine an ideal 
curve of the correlation coefficient and also of background noise. 

In a development of the lithotripter according to exemplary embodiments of the 
invention, the evaluating unit may be designed to smooth the variation of the values of the 
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correlation coefficients, for example, by averaging. The number of the correlation coefficient 
values over which the averaging operation is performed is preferably chosen to be variable. It is 
thereby prevented that, for example, an individual "outlier," which may be caused, for example, 
by a numerical artifact during calculation of a correlation coefficient, leads to the output of the 
5 error signal by the evaluating unit. Likewise, a synchronization between the ultrasonic 
transducer and the Shockwave source may be carried out to prevent electromagnetic interference 
caused by the operation of the Shockwave source and originating from one of the echoes received 
at that very moment from having a strongly interfering effect, which might lead to an outlier in 
the correlation coefficient. 

10 Furthermore, the evaluating unit may also be designed to average the minimum value 

and/or the relaxation time of the correlation coefficients over several Shockwaves. The following 
geometrical circumstances of a typical lithotripsy treatment are considered in this regard. The 
Shockwave focus has a typical extension of about 4 mm. A typical concrement, e.g. a kidney 
stone, has dimensions between 5 and 20 mm at the beginning of the treatment, and it may be 

15 moved back and forth by the mere breathing of the patient at an amplitude of about 30 mm. 
Under these conditions, a few Shockwaves would be bound to miss the concrement, because the 
emission of Shockwaves is typically carried out without regard to the respective breathing state of 
the patient. The minimum values of the correlation coefficient that are directly measured after 
such Shockwaves are particularly large, as explained above, and in an extreme case, the 

20 correlation coefficient has the constant value of 1. Accordingly, the relaxation times measured 
directly after such Shockwaves are particularly short so that, irrespective of which of the two 
parameters forms the basis for the output of an error signal, the alarm device would always be 
operated and/or the Shockwave generator would be switched off, although there is no 
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maladjustment of the patient. Instead, it is preferred that the next Shockwaves would already hit 
the concrement again because they are emitted in another breathing state of the patient. In the 
averaging of the minimum values or the relaxation time over several Shockwaves, an averaging 
operation over, for example, five Shockwaves may first be carried out in a first adjustment of the 
5 evaluating unit. But, for a further enhancement of the treatment efficiency, it is recommended 
that the number of the Shockwaves over which the averaging operation is carried out should be 
set individually at the beginning of the treatment of the patient in dependence upon his typical 
breathing behavior, the size of the concrement to be fragmented, etc. As an alternative to a 
smoothing or averaging operation, other types of signal filtering operations are also possible, 

1 0 such as a median filtering. 

In those versions of generic lithotripters that perform Doppler analyses, it has been found- 
out in clinical experiments that numerous ultrasound measurements are definitely dependent on; 
individual patients. Therefore, the evaluating unit may be designed in an expedient development 
of the lithotripter according to exemplary embodiments of the invention to standardize the 

15 minimum value and/or the relaxation time of the correlation coefficient to the reference 
minimum value or the reference relaxation time of a reference correlation coefficient curve. The 
reference correlation coefficient curve can be expediently recorded at the beginning of the 
treatment, when the patient has been adjusted in the above-described way in the lithotripter such 
that, as can be monitored with the help of the imaging locating device, a kidney stone, for 

20 instance, is located at least substantially in the Shockwave focus. In the above-mentioned display 
of the process (i.e., the continuous plotting, for example, of relaxation times) in the course of the 
Shockwave treatment, there will typically, thus, be no plotting of absolute relaxation times of the 
correlation coefficient, but standardized relaxation times will be plotted. 

14 


Atty. Docket No. 08340.Lithotripter 
(Client Docket No. PA34859US-03476/wl) 

In the procedure that has so far been described, the correlation coefficient is calculated 
according to the above formula of Eq. 1 . It has the effect that a fragmentation of the target object 
effects a sagging in the time curve of the correlation coefficient, which will relax according to a 
typical relaxation time towards the value of 1, when the system consisting of target object, 
5 fragments, and cavitation bubbles settles down again after the hit. However, the formula of Eq. 1 
will also lead to a decrease in the correlation coefficient if the concrement has, for example, 
shifted slightly due to the breathing movement of the patient between the emission of the i th 

ultrasonic pulse and the (i+k) th ultrasonic pulse. The evaluating unit operating according to the 

formula of Eq. 1 would then already calculate a decrease in the correlation coefficient, because, 

th 

10 due to the shifting of the target object, the time curve ej+k(t) of the (i+k) received ultrasonic 
echo is shifted on the time abscissa in comparison with the curve ej(t) of the i th ultrasonic echo. 
This is due to the fact that, for example, in the case of a slight movement of the concrement away 
from the ultrasonic transducer, the (i+k) th emitted ultrasonic wave has to cover a somewhat 

longer distance than the i th wave before being reflected, and as a consequence, the (i+k) th echo 

15 must travel the same additional distance. 

To eliminate such effects on the calculated correlation coefficient, it may be intended in 
some exemplary embodiments of the lithotripter that the evaluating unit is designed to determine 
a cross correlation function in time between reflected ultrasonic waves and to determine the 
maximum value of the cross correlation function in time as a correlation coefficient. In this 

20 instance the evaluating unit calculates the correlation coefficient substantially according to the 
following formula: 
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K i>k = max At 


Jei(t)e i+k (t-At)dt 


Eq.2 


The expression of Eq. 2 in brackets represents the cross correlation function with the variable At. 
Hence, in exemplary embodiments of the invention related to Eq. 2, the evaluating unit 

th 

5 shifts the measured variation of the (i+k) ultrasonic echo along the time abscissa by an amount 
At each time and it is only then that it will form the time integral, similar to the formula of Eq. 1 . 
The correlation coefficients calculated in this way for numerous values of At are temporarily 
stored, and their maximum is finally determined as the correlation coefficient that is looked 
for. This algorithm is known as a cross correlation method and is used in the medical field, for 

10 example in the determination of blood flow velocity profiles, wherein substantially that value of 
At is looked for that yields the maximum of K. Reference can be made to the publication, "Flow 
Velocity Profile via Time-Domain Correlation: Error Analysis and Computer Simulation" by 
Steven G. Foster, Paul M Embree, William D. O'Brien Jr., IEEE Transactions on Ultrasonics, 
Ferroelectrics and Frequency Control, Vol. 37, no. 2, May 1990, page 164, and to the publication, 

15 "Time Domain Formulation of Pulse-Doppler Ultrasound and Blood Velocity Estimation by 
Cross Correlation" by O. Bonnefous and P. Pesque, Ultrasonic Imaging 8, 1986, pages 73-85. 
As to possible medical information that might be contained in At, the temporary storing of the 
measuring parameters may also be implemented. 

Expediently, in the ultrasonic transducer of exemplary embodiments of the invention, the 

20 ultrasonic transmitting/receiving unit is mounted on an adjustable holder. The ultrasonic 
transducer can then be adjusted independently of other parts of the lithotripter for optimizing the 
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ultrasonic signals received and can be fixed in the optimum position. It is thereby particularly 
ensured that the ultrasonic transducer is directed to the Shockwave focus. As an alternative to the 
so-called isocentric scanner guidance, a so-called inline transducer may also be used (i.e. an 
ultrasonic transducer which is integrated into the Shockwave source). 
5 Generally, the following should be considered. With the adjustable holder, a line is just 

defined along which, in a pulse (PW) method, an emitted ultrasonic pulse propagates in the tissue 
and echoes are generated. In the case of an isocentric construction, it will then be ensured that 
the shockwave focus, i.e. the target area, is positioned along the line. With known means, such 
as a displacement sensor, it is possible to determine the distance of the transducer relative to the 

10 focus. On the basis of the known propagation time of ultrasonic pulses in the tissue, it is 
possible, in turn, to define a time window on the basis of the distance (namely, the time window 
(Ti, T 2 ) according to the formula of Eq. 1), which cuts out that part of the echo that was 
generated in the target area. 

In some exemplary embodiments of the lithotripter of the invention, the ultrasonic 

15 transmitting/receiving unit may be designed as part of an imaging ultrasonic scanner (e.g., a 
duplex scanner). In this instance, the ultrasonic transducer and parts of the electronic system of 
the ultrasonic transmitting/receiving unit can simultaneously be used for ultrasonic imaging and 
for measuring the ultrasonic echo pulses. 

Alternatively, the ultrasonic transducer may be designed as an inexpensive pin probe, 

20 which makes it flexibly usable when employed with the above-mentioned holder, for example, in 
order to ensure that it is always directed onto the Shockwave focus. Such a design is particularly 
recommended when the lithotripter, in accordance with embodiments of the invention, further 
comprises an X-ray locating device. In this instance, the imaging is carried out at the beginning 
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of the positioning of the patient and in the detail control, which takes place about every three to 
five minutes with the help of the imaging X-ray locating device, while the continuous hit control 
is carried out with the help of the ultrasonic transducer designed as the pin probe. Also in this 
instance, the correlation coefficient may be displayed on the display device that is also used for 
5 displaying the X-ray images. 

Fig. 1 is a schematic view showing major mechanical components of the lithotripter 10 
according exemplary embodiments of the invention. A patient 12 can be positioned on an 
adjustable stretcher (not shown) such that a coupling cushion 14 of a Shockwave generator 16 can 
be pressed at the desired place onto the body of the patient 12 to emit Shockwaves towards a 

10 concrement 18 to be fragmented in the body of the patient 12. In the case schematically 
illustrated in Fig. 1, the concrement 18 is a kidney stone in the kidney 20 of the patient 12. The 
patient 12 can be "adjusted" (i.e. positioned) with the help of the adjustable stretcher such that 
the focus of the Shockwaves, that are generated by a Shockwave source 22 of the Shockwave 
generator 16 and transmitted with the help of the coupling cushion 14 into the body of the patient 

15 12, is located in kidney stone 18, the focus being outlined in Fig. 1 by a cross-hair. The 
adjustment is typically carried out with the help of an imaging locating device, e.g. an X-ray 
device or an ultrasonic scanner. The components of the lithotripter 10 and an exemplary 
procedure of the positioning of a patient prior to the beginning of the treatment are generally 
known and, therefore, not be explained in further detail. 

20 Additionally, the lithotripter 10 according to exemplary embodiments of the invention can 

comprise an ultrasonic transducer 24 mounted on an adjustable holder 26. In the exemplary 
illustration of Fig. 1, the holder 26 is typically designed in the form of a hinged arm and permits 
an exact positioning of the ultrasonic transducer 25 at a desired place of the body of the patient 
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12 in such a manner that the ultrasonic transducer 24 can be directed onto the shockwave focus, 
as illustrated in Fig. 1 by dotted lines. This arrangement is called an isocentric scanner guidance. 

The ultrasonic transducer 24 transmits ultrasonic waves towards the shockwave focus in 
pulses, e.g. pulsed wave (PW), and further receives ultrasonic waves that have been reflected in 
5 the body of the patient 12, particularly from the area of the shockwave focus. As illustrated in 
Fig. 2, the received ultrasonic signals are supplied by the ultrasonic transducer 24 to a control 
unit 28, which controls the operation of the ultrasonic transducer 24 as part of a 
transmitting/receiving unit, e.g. the piezoelectric elements contained in it, and senses the 
reflected ultrasonic signals measured by the ultrasonic transducer 24. The control unit 28 also 
10 passes the signals onto downstream electronic units. For instance, the control unit 28 could pass 
on the received ultrasonic signals to an image processing module (not shown), with the help of 
which, ultrasonic images of the kidney stone 18 can be displayed on a display device 30. 

Irrespective of this feature, the control unit 28 in the lithotripter 10, according to exemplary 
1 5 embodiments of the invention, supplies the received ultrasonic signals to an evaluating unit 32, 
which is designed such that it determines a time correlation coefficient between reflected 
ultrasonic waves that are assigned to successively emitted ultrasonic pulses. These features will 
be explained with reference to Figs. 3 and 4. 

The ultrasonic transducer 24 emits ultrasonic pulses at a frequency of typically about 1 
20 kHz, which are reflected in the body of the patient 12, namely, due to the previous adjustment of 
the ultrasonic transducer 24 with the help of the adjustable holder 26 in a target area of the 
shockwave generator 16. Since the shockwave generator 16 emits about one shockwave per 
second, this means that about 1000 ultrasonic pulses are emitted into the body of the patient 12 

19 


Atty. Docket No. 08340.Lithotripter 
- '. (Client Docket No. PA34859US-03476/wl) 

between two Shockwaves and that, as a consequence, about 1000 ultrasonic echoes are also 
received. These received ultrasonic echoes are in general designated as ei, e 2 ... eiooo- Fig. 3 
shows the time-resolved characteristic of two immediately successively received ultrasonic 
echoes, ej(t) and e2(t), in a purely exemplary way, and it should be noted that the evaluating unit 
5 32 resets an internal time counter to zero again upon each emission of an ultrasonic pulse. 

In Fig. 3, the case is shown where the target area of the Shockwave generator 16 in which 
the emitted ultrasonic waves are reflected is entirely static, for example, on account of the fact 
that Shockwaves are not at all emitted at the moment, or also because of the fact that emitted 
Shockwaves miss the concrement 18. Hence, the concrement 18 remains immobile; there are 

10 neither shifts nor fragmentations of the concrement 18, and there are also no cavitation bubbles 
in the direct vicinity of the concrement 18. Accordingly, the time characteristics of the two ; 
ultrasonic echoes ei(t) and e2(t) under consideration are almost identical, so that the correlation 
calculated according to the above-indicated formula of Eq. 1 by integrating the product ei(t) • 
e2(t) over the time window of Ti to T 2 is at a maximum and will assume approximately the value 

15 of 1 upon a suitable scaling according to the above formula of Eq. 1 . 

Fig. 4 shows a case in which the ultrasonic echo ei was right away reflected on the same 
immobile target area, but the concrement 18 was hit by a Shockwave immediately thereafter. 
When the concrement 1 8 is positioned in the Shockwave focus, so that a hit is achieved upon 
emission of the Shockwave, small fragments are knocked out of the concrement 18 and the 

20 concrement 18 is simultaneously restructured internally. Furthermore, in the direct vicinity of the 
concrement 18, increased cavitation is observed in the fluid inherent to the body of patient 12. 
All of these processes have the effect that the ultrasonic waves emitted by the ultrasonic 
transducer 24 towards the concrement 18 are increasingly reflected on objects that, relative to the 
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transducer 24, are not only at another position than the one immediately before this, upon 
reflection of echo ei, but are also in movement relative to transducer 24. 

The mere fact that specific reflecting areas, for example, the concrement 1 8 itself, have 
moved based on the directly preceding reflection (e.g., away from transducer 24), leads to 
5 changed signal propagation times and thus to a shift in curve e2(t) relative to curve ei(t). 
Moreover, the fact that, due to the hit, specific reflecting areas are now no longer static, but in 
motion, effects changes in the signal curve in the reflected echo e2 (i.e. complex changes in the 
signal-time function). Additionally, the increased cavitation caused upon a hit in the direct 
vicinity of the concrement 18 can lead to an uncontrolled movement of reflecting objects (e.g., 
10 density variations or gas bubbles in the body fluid of patient 12), which can also lead to changes 
in the signal curve. 

A macroscopic movement of the concrement 1 8, which might particularly arise in cases 
where it is already fragmented at least in part, leads to both a change in signal propagation times 
and a complex change in the signal-time function when the concrement 18 performs a kind of 

15 oscillatory movement due to elastic forces of the surrounding tissue. Therefore, a look at Fig. 4 
will reveal that the correlation coefficient calculated according to the above formula of Eq. 1 will 
be smaller between the time curves ei(t) and e2(t) than in the example shown in Fig. 3, referring 
to largely identical echoes that are reflected on a static system. The evaluating unit 32 
continuously calculates the correlation coefficients between ei and e2, between e 2 and e 3 , between 

20 e3 and e4, etc., which become continuously smaller due to the above-described mechanisms in the 
case of a Shockwave hit on the concrement 18. According to experience, the fragments and the 
concrement 18 itself will settle down again after such a hit, typically after 50 to 100 ms. The 
cavitation bubbles will also have disappeared again after the end of this period, so that successive 
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echoes ej, e,+i start to resemble one another again more and more, and the correlation coefficient 
calculated by the evaluating unit 32 thus relaxes again towards its maximum value 1 (i.e. rises). 

As has already been explained above, it is preferably assumed that the characteristic of the 

f— t/T ^ 

correlation coefficient is described by a function of the form l-A K K) . However, in 
5 response to the individual situation the curve characteristic may also be similar to a Gauss curve, 
at least sectionwise. The characteristic of the correlation coefficient K to be thus expected, 
which is a function of time, is shown in Fig. 5. As described above, correlation coefficient 
values of about K = 1 can be seen at the beginning and end of the illustration, and a pronounced 
"sag" in the K curve can be seen in the middle of the illustration of Fig. 5. This sag can be 

10 described by two characteristic parameters that can be used independently and also in 
combination for identifying a hit of the Shockwave. 

On the one hand, a direct hit has the effect that the value of K falls below a predetermined 
threshold value, which is depicted in Fig. 5 by a horizontal broken line. On the other hand, it 
could generally be detected in Shockwave treatments that the relaxation time Tr (i.e., the time 

15 passing until the system of concrement 18 and of knocked-out fragments comes to rest again) is a 
measure of a hit. Hence, the width of the K curve according to the relaxation time Tr can also 
serve as a measure of a hit. Therefore, the evaluating unit 32 is designed to generate an error 
signal whenever the minimum of the correlation coefficient does not fall short of a predetermined 
first threshold value and/or whenever the relaxation time Tr of the correlation coefficient falls 

20 below a predetermined second threshold value. It can be assumed here that misses or partial hits 
of the Shockwave past the concrement 1 8 lead at best to very small sags in the K curve whereas 
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the relaxation time T R becomes longer as the concrement 18 is more precisely hit by the 
Shockwave. 

As illustrated in Fig. 2, the error signal generated by the evaluating unit 32 is supplied to 
the shockwave generator 16 and an alarm device 34. The shockwave generator 16 can thus be 
5 stopped automatically, if necessary, in case of misses. An unnecessary load on the body of the 
patient 12 by the shockwave misses is thereby avoided in a particularly reliable and efficient way 
without any intervention of the medical personnel being needed. The alarm device 34 can trigger 
an acoustic and/or optical alarm if an entered error signal indicates a miss. In such a case the 
medical personnel can manually switch off the shockwave generator 16 or take other measures, 

10 for example, newly adjusting the patient 12. 

As schematically shown in Fig. 6, a method in accordance with exemplary embodiments 
of the invention includes three sections after the lithotripter 10 has been started, namely a first 
section S10 with steps for pre-adjusting the lithotripter before the therapy proper, a subsequent 
section S20 in which, during therapy, the relevant relaxation time Tr is determined in various 

15 steps of the method, and a further section S30 in which the corresponding measures are taken in 
several method steps as a function of the relevant relaxation time T R determined in section S20. 
It should be noted that for the sake of simplicity, only the sensed relaxation time T R is used in the 
following as a basis for the hit control, but it is of course also possible that the minimum of the 
correlation coefficient K according to Fig. 5 or a combination of both parameters is used by the 

20 evaluating unit 32 for generating an error signal. 

Exemplary steps of the first section S10 can be explained with reference to Fig. 7. First, 
the patient is adjusted in step Sll with the help of the adjustable stretcher of the lithotripter 10. 
At the end of this adjustment, which is monitored by an imaging locating device, the patient 12 is 
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positioned such that the concrement 18 is within the Shockwave focus. Subsequently, the 
emission of ultrasonic waves by the ultrasonic transducer 24 and the emission of a shockwave for 
the first time with the help of shockwave generator 16 are started in step SI 2. 

Subsequently, in step SI 3, the evaluating unit 32 determines the characteristic of the 
5 correlation coefficient with the help of the reflected ultrasonic echo in the above-described way, 
the correlation coefficient being then displayed in step S14 on the display device 30 according to 
Fig. 5. A fit curve, which has been adapted (i.e., "fitted") by the evaluating unit 32 to the sag in 
the K curve can already be displayed on the display device 30 in this step. This characteristic of 
the correlation coefficient is inspected in step S15 by the medical personnel as to whether it 

1 0 indicates a hit or a miss. 

Typically, this first emission of a shockwave in step S12 yields a hit, because after step 
SI 1, the concrement 18 is positioned within the shockwave focus, and because prior to the 
emission of this one "test shockwave," one can ask the patient 1 2 to stop breathing for a short 
period of time, so that a miss due to a breathing movement of the concrement 1 8 is not imminent. 

15 In this inspection of correlation coefficient characteristics as to whether there is a hit or a miss, it 
is, among other things, important that the respective signal noise is determined, such as 
variations that arise upon reflection of the ultrasonic wave on a resting patient 12 without the 
emission of Shockwaves. If a hit is not detected in step SI 5, step SI 1 is resumed. By contrast, if 
a hit is detected, the measured correlation coefficient curve is stored in step SI 6 as a reference 

20 correlation coefficient curve in memory 32a of the evaluating unit 32. 

With the help of the reference correlation coefficient curve in memory 32a, the above- 
discussed threshold values are then defined in step SI 7. In case these are exceeded or not 
reached, this will have the effect that the evaluating unit 32 supplies an error signal. 
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Furthermore, in step SI 7, the relaxation time T R , which is determined in the reference correlation 
coefficient curve, is stored in memory 32a as a reference relaxation time. An initialization (i.e. a 
resetting of a count parameter Z), whose significance will become clear in the following from 
Fig. 8, is then performed in step S 1 8. 
5 Fig. 8 shows exemplary steps within the second section S20 according to Fig. 6 which 

may serve to determine the relevant relaxation time T R . First, the count parameter Z is 
incremented in step S21 (e.g., increased by 1). A Shockwave is subsequently emitted in step S22 
by the Shockwave generator 26. The characteristic of the correlation coefficient is determined in 
step S23 and subsequently displayed in step S24 on the display device 30, just like an optionally 

1 0 adapted fit curve. 

In the subsequent step S25, the relaxation time Tr of the actual correlation coefficient, 
curve is determined with the help of the fitted Gauss curve (other fit curves or evaluation 
methods are possible, such as those fundamentally known in the field of data evaluation). In step 
S26, the relaxation time T R is converted by means of a division by the reference relaxation time 

15 determined in step SI 7, according to Fig. 7, into a standardized relaxation time Tr, which will 
then be stored in step S27 in the evaluating unit 42, for example, in the memory 32a or in a 
separate memory for receiving measured values. 

It is then checked, in step S28, whether the count parameter Z has reached a 
predetermined value, namely the value Z = 5 in the example shown in Fig. 8. If this is not the 

20 case, the program executed in the evaluating unit 32 returns to step S21, in which the count 
parameter Z is incremented and a further correlation coefficient characteristic is then measured. 
However, when the examination in step S28 reveals that the count parameter Z has reached the 
predetermined value, for example Z = 5 (which means that the relaxation times of 5 preceding 
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correlation coefficient characteristics are temporarily stored in the evaluating unit 32), the 
program continues and goes to a step S29 in which an averaging operation over the five 
temporarily stored relaxation times is carried out. The mean value of the relaxation times of the 
last five correlation coefficient curves that are calculated in step S29 can also be stored 
5 temporarily in the evaluating unit 32 in one of the memories. 

In step S29a, the mean value that is the last calculated one is added to a continuous 
display, which continuously displays the relaxation times measured in the course of the whole 
treatment to trace the development of the relaxation times. Expediently, the continuous display 
can be shown on the display device 30, but this may also be a separate screen, etc. An example 

10 of such a continuous display is shown in Fig. 10 and shall be explained further below. 

Fig. 9 shows exemplary steps of the third section S30 according to Fig. 6 and the 
measures that can be taken within the scope of the method of the invention during a lithotripsy 
treatment with the lithotripter 10 of the invention as a function of the mean value calculated in 
step S29. In a step S31, it is checked through the program executed in the evaluating unit 32 

15 whether the mean value calculated in step S29 is larger than a predetermined threshold value. As 
explained above, large relaxation time values indicate a hit, whereas small relaxation time values 
are typical of misses. Consequently, a positive result of the examination in step S31 means that 
the concrement 18 is apparently still within the Shockwave focus and hits are thus made. As a 
result in this case, the program executed in evaluating unit 32 goes back to S21 via step S31, in 

a 

20 which the count parameter Z is reset to 0 (i.e. the program starts a new measurement of five 
standardized relaxation times with a subsequent averaging operation according to Fig. 8). 
However, when the examination in step S3 1 reveals a negative result, for example, the averaged 
relaxation time T R of the last five correlation coefficient characteristics is too short, (which 
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indicates possible misses), the evaluating unit 32 outputs an error signal and the program 
branches off to a step S33 in which the alarm device 34 is triggered, and subsequently to a step 
S34 in which the Shockwave generator 16 is stopped to avoid unnecessary load on the body of 
the patient 1 2 by misses. 

5 Fig. 10 is a schematic illustration showing an exemplary variation control display on the 

display device 30 during a Shockwave therapy. Fig. 10 reveals a coordinate system, the abscissa 
of which indicates the number of emitted Shockwaves, and the ordinate of which indicates the 
absolute or standardized relaxation time T R determined subsequent to the respective Shockwave. 
As depicted, from the start of the Shockwave therapy until shortly before the end thereof, there is 

10 a characteristic corresponding to a slow rise in relaxation time Tr. Individual "sags" in this 
characteristic curve indicate misses that have been caused by a movement of the patient 12 or of 
the concrement 18 within the patient 12 and that have been corrected immediately by newly 
positioning the patient 12. In Fig. 9, a horizontal broken line represents the amount of the 
threshold value which will trigger the alarm device 34 if not reached. As can be seen, the local 

15 sag has led approximately in the middle of the characteristic curve to a falling short of the 
threshold value and thus to an alarm. 

It is noted that Figs. 6-10 show exemplary steps of methods in accordance with exemplary 
embodiments of the invention and that numerous additional preceding, succeeding, and/or 
intermediate steps are possible, including those that are basically known in the field of lithotripsy 

20 methods. It also noted that, as mentioned above, the minimum value of the correlation 
coefficient (e.g., in Fig. 5) can also be used instead of or in addition to the relaxation time T R to 
detect a miss. Since in the illustration of Fig. 5 hits lead to deep and simultaneously broad sags, 
the area contained in the "sag" of the K curve can also be determined numerically and used as a 
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measure of a hit or a miss. Furthermore, it is noted that the number Z = 5, as indicated in step 
S28, is a mere exemplary value. To be more specific, this number can also be entered 
individually by medical personnel within the pre-adjustments according to section S10 in Fig. 6 
(e.g. in step SI 7 in Fig. 7). Hence, in a patient 12 with a rather small concrement 18 (which can 
5 thus be missed by the shockwave more easily) and heavy breathing at the same time, with 
correspondingly strong deflections of the concrement 18, the averaging operation could be 
performed over a larger number of Shockwaves and, thus, relaxation times than in another patient 
12 with a rather flat breathing and a larger concrement 18 (which is thus hit more often). It is 
also possible to work without any averaging operation (e.g., Z = 1). 

10 It is further noted that the reference relaxation time (step SI 7) and the later 

standardization of measured relaxation times, with the help of the reference relaxation time (step 
S26), need not be determined so that in the variation control during shockwave therapy absolute 
values of the measured relaxation times are displayed on the display device 20 and evaluated 
accordingly. As for the lithotripter 10, which is outlined schematically in Fig. 1, it should be 

15 understood that the device 10 may comprise numerous further components, including those that 
are known in the art, such as an X-ray locating device or an imaging ultrasonic scanner. 

While the invention has been described with respect to the foregoing exemplary 
embodiments, it will be apparent to those skilled in the art that various modifications, variations 
and improvements of the invention may be made in light of the above teachings and within the 

20 purview of the appended claims without departing from the spirit and intended scope of the 
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invention. In regard to the foregoing description of the exemplary embodiments of the invention, 
areas which are known to those of ordinary skill in the art may not have been described in detail, 
in order to facilitate a clear and concise description of the invention. Accordingly, it should be 
understood that the invention is not to be limited by the specific exemplary embodiments, but 
only by the scope of the appended claims. 
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